QTS data: Spiritual curiosity relationships

Last updated 2018-04-12

This is a first-pass exploration of some of the relationships among QTS data for the “Spiritual curiosity” interview.

Time praying & spiritual experiences

Requested by Josh

Time praying as indexed by prayfreqmin

How many minutes each day do you do that on average throughout the week? [use real world examples, ask alone? With others? Probe away– Answer with interviewer judgement for answers]

Spiritual experiences indexed by godvoxaloud

Some people say that they have heard God* speak out loud to them. Has this ever happened to you?

Ignoring unknown aesthetics: fill

Ignoring unknown aesthetics: fill

Call:
glm(formula = godvoxaloud_num ~ prayfreqmin + country, family = "binomial", 
    data = d1)

Deviance Residuals: 
   Min      1Q  Median      3Q     Max  
-1.703  -1.027  -0.455   1.032   2.154  

Coefficients:
              Estimate Std. Error z value Pr(>|z|)    
(Intercept)   -0.47926    0.13784  -3.477 0.000507 ***
prayfreqmin.L  0.73747    0.24453   3.016 0.002562 ** 
prayfreqmin.Q  0.54372    0.27178   2.001 0.045436 *  
prayfreqmin.C -0.35794    0.29205  -1.226 0.220347    
countryGH      0.97536    0.21961   4.441 8.94e-06 ***
countryTH      0.25817    0.26628   0.970 0.332274    
countryCH      0.30449    0.24600   1.238 0.215809    
countryVT      0.05578    0.25920   0.215 0.829609    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 446.17  on 328  degrees of freedom
Residual deviance: 385.93  on 321  degrees of freedom
AIC: 401.93

Number of Fisher Scoring iterations: 4

Spiritual experiences indexed by godviavisions

_Some people say that they have had a vision from God*—they have a picture, but it is like they see it with their eyes. Has anything like that happened to you?_

Ignoring unknown aesthetics: fill

Ignoring unknown aesthetics: fill

Call:
glm(formula = godviavisions_num ~ prayfreqmin + country, family = "binomial", 
    data = d2)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.8640  -1.0136  -0.5422   1.1679   2.0691  

Coefficients:
              Estimate Std. Error z value Pr(>|z|)    
(Intercept)    -0.6210     0.1361  -4.562 5.06e-06 ***
prayfreqmin.L   0.9386     0.2470   3.800 0.000145 ***
prayfreqmin.Q   0.8468     0.2739   3.092 0.001989 ** 
prayfreqmin.C   0.3692     0.2926   1.262 0.207031    
countryGH       1.0290     0.2239   4.597 4.29e-06 ***
countryTH       0.6839     0.2719   2.515 0.011894 *  
countryCH      -0.9077     0.2933  -3.095 0.001968 ** 
countryVT       0.1282     0.2570   0.499 0.617930    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 443.63  on 334  degrees of freedom
Residual deviance: 386.58  on 327  degrees of freedom
AIC: 402.58

Number of Fisher Scoring iterations: 4

Spiritual experiences indexed by godviabodyexperiences

Some people have particular experiences in your body that they associate with God* or spirit. Does that happen for you? [examples: warm hands, goosebumps, fire in the belly]

Ignoring unknown aesthetics: fill

Ignoring unknown aesthetics: fill

Call:
glm(formula = godviabodyexperiences_num ~ prayfreqmin + country, 
    family = "binomial", data = d3)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-2.0197  -1.1520   0.6754   0.8918   1.3475  

Coefficients:
              Estimate Std. Error z value Pr(>|z|)    
(Intercept)     0.7061     0.1332   5.302 1.14e-07 ***
prayfreqmin.L   0.5249     0.2378   2.207  0.02731 *  
prayfreqmin.Q   0.5157     0.2622   1.966  0.04926 *  
prayfreqmin.C   0.2188     0.2779   0.788  0.43096    
countryGH       0.5353     0.2375   2.254  0.02418 *  
countryTH       0.8843     0.3073   2.877  0.00401 ** 
countryCH      -0.4790     0.2488  -1.926  0.05415 .  
countryVT      -0.1303     0.2500  -0.521  0.60225    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 432.22  on 334  degrees of freedom
Residual deviance: 402.36  on 327  degrees of freedom
AIC: 418.36

Number of Fisher Scoring iterations: 4

Aggregate spiritual experiences score

Linear mixed model fit by REML ['lmerMod']
Formula: response ~ prayfreqmin + country + (1 | subject_name) + (1 |  
    question)
   Data: d_agg1

REML criterion at convergence: 9516

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.6638 -0.7863 -0.1873  0.8822  2.5579 

Random effects:
 Groups       Name        Variance Std.Dev.
 subject_name (Intercept) 0.02090  0.1446  
 question     (Intercept) 0.03517  0.1875  
 Residual                 0.17972  0.4239  
Number of obs: 7961, groups:  subject_name, 344; question, 24

Fixed effects:
               Estimate Std. Error t value
(Intercept)    0.459403   0.039550  11.616
prayfreqmin.L  0.117263   0.017788   6.592
prayfreqmin.Q  0.033800   0.019728   1.713
prayfreqmin.C  0.048443   0.021485   2.255
countryGH      0.024153   0.016975   1.423
countryTH      0.001581   0.020844   0.076
countryCH     -0.086277   0.019517  -4.421
countryVT      0.155426   0.019417   8.005

Correlation of Fixed Effects:
            (Intr) pryf.L pryf.Q pryf.C cntrGH cntrTH cntrCH
prayfrqmn.L  0.035                                          
prayfrqmn.Q -0.066  0.174                                   
prayfrqmn.C  0.024 -0.164 -0.026                            
countryGH   -0.045 -0.099  0.114  0.118                     
countryTH    0.034  0.240  0.040 -0.020 -0.259              
countryCH    0.014 -0.076 -0.095 -0.034 -0.227 -0.305       
countryVT    0.004 -0.045  0.058 -0.022 -0.209 -0.299 -0.278
Linear mixed model fit by REML ['lmerMod']
Formula: response ~ prayfreqmin * country + (1 | subject_name) + (1 |  
    question)
   Data: d_agg1

REML criterion at convergence: 9563.5

Scaled residuals: 
    Min      1Q  Median      3Q     Max 
-2.6672 -0.7887 -0.1796  0.8787  2.5708 

Random effects:
 Groups       Name        Variance Std.Dev.
 subject_name (Intercept) 0.02118  0.1455  
 question     (Intercept) 0.03517  0.1875  
 Residual                 0.17973  0.4239  
Number of obs: 7961, groups:  subject_name, 344; question, 24

Fixed effects:
                         Estimate Std. Error t value
(Intercept)              0.458378   0.039986  11.463
prayfreqmin.L            0.108947   0.022863   4.765
prayfreqmin.Q            0.026234   0.023102   1.136
prayfreqmin.C            0.042343   0.023338   1.814
countryGH                0.027652   0.018607   1.486
countryTH               -0.012737   0.030968  -0.411
countryCH               -0.075852   0.022478  -3.374
countryVT                0.153817   0.020838   7.381
prayfreqmin.L:countryGH  0.006461   0.035854   0.180
prayfreqmin.Q:countryGH  0.003627   0.037213   0.097
prayfreqmin.C:countryGH  0.023226   0.038524   0.603
prayfreqmin.L:countryTH -0.039139   0.069967  -0.559
prayfreqmin.Q:countryTH -0.016074   0.061933  -0.260
prayfreqmin.C:countryTH  0.009706   0.052689   0.184
prayfreqmin.L:countryCH -0.029814   0.038224  -0.780
prayfreqmin.Q:countryCH -0.039246   0.044956  -0.873
prayfreqmin.C:countryCH  0.018775   0.050805   0.370
prayfreqmin.L:countryVT -0.012807   0.039717  -0.322
prayfreqmin.Q:countryVT  0.036971   0.041676   0.887
prayfreqmin.C:countryVT  0.035715   0.043548   0.820

Correlation matrix not shown by default, as p = 20 > 12.
Use print(x, correlation=TRUE)  or
     vcov(x)     if you need it

refitting model(s) with ML (instead of REML)

Data: d_agg1
Models:
r_agg1: response ~ prayfreqmin + country + (1 | subject_name) + (1 | 
r_agg1:     question)
r_agg2: response ~ prayfreqmin * country + (1 | subject_name) + (1 | 
r_agg2:     question)
       Df    AIC    BIC  logLik deviance Chisq Chi Df Pr(>Chisq)
r_agg1 11 9489.8 9566.6 -4733.9   9467.8                        
r_agg2 23 9504.7 9665.3 -4729.3   9458.7 9.134     12     0.6914

Ignoring unknown aesthetics: y

Looking just at ‘apples’ (urban charismatics)

Looking at all charismatics

Doubt & external sensory experiences

Requested by Emily, John

“Doubt” indxed by selfunsuregodreal

Has there been a time when you yourself wondered whether God* was real?

Spiritual experiences indexed by godvoxaloud

Some people say that they have heard God* speak out loud to them. Has this ever happened to you?

Ignoring unknown aesthetics: fill

Ignoring unknown aesthetics: fill

Call:
glm(formula = godvoxaloud_num ~ selfunsuregodreal + country, 
    family = "binomial", data = d4)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.4270  -1.0266  -0.5155   0.9986   2.0415  

Coefficients:
                   Estimate Std. Error z value Pr(>|z|)    
(Intercept)        -0.47400    0.12983  -3.651 0.000261 ***
selfunsuregodrealY -0.06683    0.13151  -0.508 0.611315    
countryGH           0.97704    0.21634   4.516 6.29e-06 ***
countryTH           0.01075    0.26818   0.040 0.968028    
countryCH           0.38097    0.24572   1.550 0.121032    
countryVT           0.04141    0.26124   0.159 0.874046    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 426.65  on 313  degrees of freedom
Residual deviance: 381.81  on 308  degrees of freedom
AIC: 393.81

Number of Fisher Scoring iterations: 4

Spiritual experiences indexed by godviavisions

_Some people say that they have had a vision from God*—they have a picture, but it is like they see it with their eyes. Has anything like that happened to you?_

Ignoring unknown aesthetics: fill

Ignoring unknown aesthetics: fill

Call:
glm(formula = godviavisions_num ~ selfunsuregodreal + country, 
    family = "binomial", data = d5)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.3530  -0.9346  -0.6884   1.0629   1.8214  

Coefficients:
                   Estimate Std. Error z value Pr(>|z|)    
(Intercept)        -0.61222    0.12652  -4.839 1.30e-06 ***
selfunsuregodrealY  0.06415    0.13016   0.493   0.6221    
countryGH           0.95182    0.21279   4.473 7.71e-06 ***
countryTH           0.33395    0.26614   1.255   0.2096    
countryCH          -0.58890    0.27649  -2.130   0.0332 *  
countryVT           0.07429    0.25358   0.293   0.7696    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 422.46  on 318  degrees of freedom
Residual deviance: 391.05  on 313  degrees of freedom
AIC: 403.05

Number of Fisher Scoring iterations: 4

Spiritual experiences indexed by godviabodyexperiences

Some people have particular experiences in your body that they associate with God* or spirit. Does that happen for you? [examples: warm hands, goosebumps, fire in the belly]

Ignoring unknown aesthetics: fill

Ignoring unknown aesthetics: fill

Call:
glm(formula = godviabodyexperiences_num ~ selfunsuregodreal + 
    country, family = "binomial", data = d6)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.8606  -1.2321   0.6274   0.9309   1.2922  

Coefficients:
                   Estimate Std. Error z value Pr(>|z|)    
(Intercept)          0.7056     0.1277   5.523 3.33e-08 ***
selfunsuregodrealY   0.2421     0.1317   1.839   0.0659 .  
countryGH            0.5778     0.2363   2.446   0.0145 *  
countryTH            0.5884     0.3089   1.905   0.0568 .  
countryCH           -0.3357     0.2470  -1.359   0.1742    
countryVT           -0.1013     0.2519  -0.402   0.6876    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 411.83  on 319  degrees of freedom
Residual deviance: 390.07  on 314  degrees of freedom
AIC: 402.07

Number of Fisher Scoring iterations: 4

“Doubt” indxed by morequesmoreanswr

Do you think that the more spiritually mature you become, you will discover more questions or more answers?

Spiritual experiences indexed by godvoxaloud

Some people say that they have heard God* speak out loud to them. Has this ever happened to you?

Ignoring unknown aesthetics: fill

Ignoring unknown aesthetics: fill

Call:
glm(formula = godvoxaloud_num ~ morequesmoreanswr + country, 
    family = "binomial", data = d7)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.4044  -0.9984  -0.4288   0.9900   2.2055  

Coefficients:
                   Estimate Std. Error z value Pr(>|z|)    
(Intercept)        -0.65574    0.16797  -3.904 9.46e-05 ***
morequesmoreanswrQ -0.03057    0.15819  -0.193    0.847    
countryGH           1.14447    0.23706   4.828 1.38e-06 ***
countryTH           0.18839    0.31250   0.603    0.547    
countryCH           0.28946    0.33757   0.857    0.391    
countryVT           0.03154    0.32558   0.097    0.923    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 325.90  on 241  degrees of freedom
Residual deviance: 280.75  on 236  degrees of freedom
AIC: 292.75

Number of Fisher Scoring iterations: 4

Spiritual experiences indexed by godviavisions

_Some people say that they have had a vision from God*—they have a picture, but it is like they see it with their eyes. Has anything like that happened to you?_

Ignoring unknown aesthetics: fill

Ignoring unknown aesthetics: fill

Call:
glm(formula = godviavisions_num ~ morequesmoreanswr + country, 
    family = "binomial", data = d8)

Deviance Residuals: 
   Min      1Q  Median      3Q     Max  
-1.388  -1.039  -0.616   1.127   2.026  

Coefficients:
                   Estimate Std. Error z value Pr(>|z|)    
(Intercept)         -0.6647     0.1687  -3.939 8.18e-05 ***
morequesmoreanswrQ   0.1811     0.1525   1.188   0.2349    
countryGH            0.9661     0.2353   4.106 4.03e-05 ***
countryTH            0.5124     0.3171   1.616   0.1061    
countryCH           -1.0697     0.4511  -2.372   0.0177 *  
countryVT            0.3111     0.3162   0.984   0.3251    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 325.22  on 242  degrees of freedom
Residual deviance: 295.17  on 237  degrees of freedom
AIC: 307.17

Number of Fisher Scoring iterations: 4

Spiritual experiences indexed by godviabodyexperiences

Some people have particular experiences in your body that they associate with God* or spirit. Does that happen for you? [examples: warm hands, goosebumps, fire in the belly]

Ignoring unknown aesthetics: fill

Ignoring unknown aesthetics: fill

Call:
glm(formula = godviabodyexperiences_num ~ morequesmoreanswr + 
    country, family = "binomial", data = d9)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-1.9152  -1.2161   0.7310   0.9641   1.1674  

Coefficients:
                   Estimate Std. Error z value Pr(>|z|)    
(Intercept)         0.71079    0.15975   4.449 8.62e-06 ***
morequesmoreanswrQ  0.03348    0.15584   0.215   0.8299    
countryGH           0.43886    0.24850   1.766   0.0774 .  
countryTH           0.91566    0.38084   2.404   0.0162 *  
countryCH          -0.54853    0.33456  -1.640   0.1011    
countryVT          -0.15232    0.30627  -0.497   0.6190    
---
Signif. codes:  0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1

(Dispersion parameter for binomial family taken to be 1)

    Null deviance: 310.96  on 244  degrees of freedom
Residual deviance: 294.81  on 239  degrees of freedom
AIC: 306.81

Number of Fisher Scoring iterations: 4

To do

Here are other things on our to-do list:

• Experiences by sense (Josh)
• Extremity of experiences (Josh)
• …